X 



Measurement of branching ratio and B® lifetime in the decay B° s — > J/ipf b(980) at CDF 

T. Aaltonen, 21 B. Alvarez Gonzalez 1 ", 9 S. Amerio, 41 D. Amidei, 32 A. Anastassov, 36 A. Annovi, 17 J. Antos, 12 
G. Apollinari, 15 J. A. Appel, 15 A. Apresyan, 46 T. Arisawa, 56 A. Artikov, 13 J. Asaadi, 51 W. Ashmanskas, 15 
B. Auerbach, 59 A. Aurisano, 51 F. Azfar, 40 W. Badgett, 15 A. Barbaro-Galtieri, 26 V.E. Barnes, 46 B.A. Barnett, 23 
P. Barria dd , 44 P. Bartos, 12 M. Bauce bh , 41 G. Bauer, 30 F. Bedeschi, 44 D. Beecher, 28 S. Behari, 23 G. Bellettini cc , 44 
J. Bellinger, 58 D. Benjamin, 14 A. Beretvas, 15 A. Bhatti, 48 M. Binkley*, 15 D. Bisello 66 , 41 I. Bizjak^, 28 K.R. Bland, 5 
B. Blumcnfeld. 23 A. Bocci 14 A. Bodek. 47 D. Bortoletto. 46 .7. Boudreani. 45 A. Boveia. 11 B. Brau a 15 



nger,°° D. Benjamin, 1 * A. Beretvas, 10 A. Bhatti,"" M. Binkley*, 10 D. Bisello°V 1. Bizjak" 
B. Blumenfeld, 23 A. Bocci, 14 A. Bodek, 47 D. Bortoletto, 46 J. Boudreau, 45 A. Boveia, 11 B. Brau a , 15 
L. Brigliadori aa , 6 A. Brisuda, 12 C. Bromberg, 33 E. Brucken, 21 M. Bucciantonio cc , 44 J. Budagov, 13 
H.S. Budd, 47 S. Budd, 22 K. Burkett, 15 G. Busetto 66 , 41 P. Bussey, 19 A. Buzatu, 31 C. Calancha, 29 S. Camarda, 

M. Campanelli, 33 M. Campbell, 32 F. Canelli 11 , 15 B. Carls, 22 D. Carlsmith, 58 R. Carosi, 44 S. CarrilloY 6 
3. Carron, 15 B. Casal, 9 M. Casarsa, 15 A. Castro QQ , 6 P. Catastini, 20 D. Cauz, 52 V. Cavaliere, 22 M. Cavalli-Sforz 

4 florri/ 26 T. Porritn? 2 » V C P1,pti 1 M nViorfr.L- 7 CI Ptiiarolli 44 (~< rMxr-hUvo 15 p aioKana « K r.Yir* '■ 



O 

(N ; c 



JVL. Campanelli, 00 M. Campbell, oz F. Canelli 11 , 10 B. Carls, zz D. Carlsmith, 00 R. Carosi,^ S. CarrilloV 
!. Carron, 15 B. Casal, 9 M. Casarsa, 15 A. Castro QQ , 6 P. Catastini, 20 D. Cauz, 52 V. Cavaliere, 22 M. Cavalli-Sforza, 4 
A. Cerri^, 26 L. Cerrito 9 , 28 Y.C. Chen, 1 M. Chertok, 7 G. Chiarelli, 44 G. Chlachidze, 15 F. Chlebana, 15 K. Cho, 25 

D. Chokhcli, 13 J.P. Chou, 20 W.H. Chung, 58 Y.S. Chung, 47 C.I. Ciobanu, 42 M.A. Ciocci drf , 44 A. Clark, 18 
C. Clarke, 57 G. Compostella bh , 41 M.E. Convery, 15 J. Conway, 7 M.Corbo, 42 M. Cordelli, 17 C.A. Cox, 7 D.J. Cox, 7 

F. Crescioli cc , 44 C. Cuenca Almenar, 59 J. Cuevas u ', 9 R. Culbertson, 15 D. Dagenhart, 15 N. d'Ascenzo", 42 
M. Datta, 15 P. de Barbara, 47 S. De Cecco, 49 G. De Lorenzo, 4 M. Dell'Orso cc , 44 C. Deluca, 4 L. Demortier, 48 

J. Deng c , 14 M. Deninno, 6 F. Devoto, 21 M. d'Errico h& , 41 A. Di Canto cc , 44 B. Di Ruzza, 44 J.R. Dittmann, 5 
M. D'Onofrio, 27 S. Donati cc , 44 P. Dong, 15 M. Dorigo, 52 T. Dorigo, 41 K. Ebina, 56 A. Elagin, 51 A. Eppig, 32 
R. Erbacher, 7 D. Errede, 22 S. Errede, 22 N. Ershaidat 2 , 42 R. Eusebi, 51 H.C. Fang, 26 S. Farrington, 40 M. Feindt, 24 
J.P. Fernandez, 29 C. Ferrazza ee , 44 R. Field, 16 G. Flanagan 3 , 46 R. Forrest, 7 M.J. Frank, 5 M. Franklin, 20 
J.C. Freeman, 15 Y. Funakoshi, 56 I. Furic, 16 M. Gallinaro, 48 J. Galyardt, 10 J.E. Garcia, 18 A.F. Garfinkel, 46 
P. Garosi'", 44 H. Gerberich, 22 E. Gerchtein, 15 S. Giagu // , 49 V. Giakoumopoulou, 3 P. Giannetti, 44 K. Gibson, 45 

CM. Ginsburg, 15 N. Giokaris, 3 P. Giromini, 17 M. Giunta, 44 G. Giurgiu, 23 V. Glagolev, 13 D. Glenzinski, 15 
M. Gold, 35 D. Goldin, 51 N. Goldschmidt, 16 A. Golossanov, 15 G. Gomez, 9 G. Gomez-Ceballos, 30 M. Goncharov, 30 
O. Gonzalez, 29 I. Gorelov, 35 A.T. Goshaw, 14 K. Goulianos, 48 S. Grinstein, 4 C. Grosso-Pilcher, 11 R.C. Group 55 , 15 
J. Guimaraes da Costa, 20 Z. Gunay-Unalan, 33 C. Haber, 26 S.R. Hahn, 15 E. Halkiadakis, 50 A. Hamaguchi, 39 
JY. Han, 47 F. Happacher, 17 K. Hara, 53 D. Hare, 50 M. Hare, 54 R.F. Harr, 57 K. Hatakeyama, 5 C. Hays, 40 
M. Heck, 24 J. Heinrich, 43 M. Herndon, 58 S. Hewamanage, 5 D. Hidas, 50 A. Hocker, 15 W. Hopkins 9 , 15 D. Horn, 24 
S. Hou, 1 R.E. Hughes, 37 M. Hurwitz, 11 M. Huschle, 24 U. Husemann, 59 N. Hussain, 31 M. Hussein, 33 J. Huston, 33 

G. Introzzi, 44 M. Iori^, 49 A. Ivanov , 7 E. James, 15 D. Jang, 10 B. Jayatilaka, 14 E.J. Jcon, 25 M.K. Jha, 6 
S. Jindariani, 15 W. Johnson, 7 M. Jones, 46 K.K. Joo, 25 SY. Jun, 10 T.R. Junk, 15 T. Kamon, 51 P.E. Karchin, 57 

A. Kasmi, 5 Y. Kato™, 39 W. Ketchum, 11 J. Kcung, 43 V. Khotilovich, 51 B. Kilminster, 15 D.H. Kim, 25 H.S. Kim, 25 

H. W. Kim, 25 J.E. Kim, 25 M.J. Kim, 17 S.B. Kim, 25 S.H. Kim, 53 Y.K. Kim, 11 N. Kimura, 56 M. Kirby, 15 
S. Klimenko, 16 K. Kondo, 56 D.J. Kong, 25 J. Konigsberg, 16 A.V. Kotwal, 14 M. Kreps", 24 J. Kroll, 43 D. Krop, 11 

N. Krumnack', 5 M. Kruse, 14 V. Krutelyov d , 51 T. Kuhr, 24 M. Kurata, 53 S. Kwang, 11 A.T. Laasanen, 46 S. Lami, 44 

S. Lammel, 15 M. Lancaster, 28 R.L. Lander, 7 K. Lannon 11 , 37 A. Lath, 50 G. Latino cc , 44 T. LeCompte, 2 E. Lee, 51 
H.S. Lee, 11 J.S. Lee, 25 S.W. Lee*, 51 S. Leo cc , 44 S. Leone, 44 J.D. Lewis, 15 A. Limosani r , 14 C.-J. Lin, 26 J. Linacre, 40 

M. Lindgren, 15 E. Lipeles, 43 A. Lister, 18 D.O. Litvintsev, 15 C. Liu, 45 Q. Liu, 46 T. Liu, 15 S. Lockwitz, 59 
A. Loginov, 59 D. Lucchesi &fc , 41 J. Lueck, 24 P. Lujan, 26 P. Lukens, 15 G. Lungu, 48 J. Lys, 26 R. Lysak, 12 R. Madrak, 15 
K. Maeshima, 15 K. Makhoul, 30 S. Malik, 48 G. Manca 6 , 27 A. Manousakis-Katsikakis, 3 F. Margaroli, 46 C. Marino, 24 

M. Martinez, 4 R. Martmez-Ballarm, 29 P. Mastrandrea, 49 M.E. Mattson, 57 P. Mazzanti, 6 K.S. McFarland, 47 
P. Mclntyre, 51 R. McNulty*, 27 A. Mehta, 27 P. Mehtala, 21 A. Menzione, 44 C. Mesropian, 48 T. Miao, 15 D. Mietlicki, 32 

A. Mitra, 1 H. Miyake, 53 S. Moed, 20 N. Moggi, 6 M.N. Mondragon fc , 15 C.S. Moon, 25 R. Moore, 15 M.J. Morello, 15 
J. Morlock, 24 P. Movilla Fernandez, 15 A. Mukherjee, 15 Th. Muller, 24 P. Murat, 15 M. Mussini aa , 6 J. Nachtman™, 15 
Y. Nagai, 53 J. Naganoma, 56 I. Nakano, 38 A. Napier, 54 J. Nett, 51 C. Neu, 60 M.S. Neubauer, 22 J. Nielsen 6 , 26 
L. Nodulman, 2 O. Norniella, 22 E. Nurse, 28 L. Oakes, 40 S.H. Oh, 14 Y.D. Oh, 25 I. Oksuzian, 60 T. Okusawa, 39 
R. Orava, 21 L. Ortolan, 4 S. Pagan Griso hb , 41 C. Pagliarone, 52 E. Palencia^, 9 V. Papadimitriou, 15 A. A. Paramonov, 2 
J. Patrick, 15 G. Pauletta 99 , 52 M. Paulini, 10 C. Paus, 30 D.E. Pellett, 7 A. Penzo, 52 T.J. Phillips, 14 G. Piacentino, 44 
E. Pianori, 43 J. Pilot, 37 K. Pitts, 22 C. Plager, 8 L. Pondrom, 58 K. Potamianos, 46 O. Poukhov*, 13 F. Prokoshin 9 , 13 

A. Pronko, 15 F. Ptohos' 1 , 17 E. Pueschel, 10 G. Punzi cc , 44 J. Pursley, 58 A. Rahaman, 45 V. Ramakrishnan, 58 
N. Ranjan, 46 I. Redondo, 29 P. Renton, 40 M. Rescigno, 49 T. Riddick, 28 F. Rimondi QQ , 6 L. Ristori 44 , 15 A. Robson, 19 
T. Rodrigo, 9 T. Rodriguez, 43 E. Rogers, 22 S. Rolli, 54 R. Roser, 15 M. Rossi, 52 F. Rubbo, 15 F. Ruffini dti , 44 
A. Ruiz, 9 J. Russ, 10 V. Rusu, 15 A. Safonov, 51 W.K. Sakumoto, 47 Y. Sakurai, 56 L. Santi 99 , 52 L. Sartori, 44 



2 



K. Sato, 53 V. Saveliev 11 , 42 A. Savoy-Navarro, 42 P. Schlabach, 15 A. Schmidt, 24 E.E. Schmidt, 15 M.P. Schmidt*, 59 
M. Schmitt, 36 T. Schwarz, 7 L. Scodellaro, 9 A. Scribano'", 44 F. Scuri, 44 A. Scdov, 46 S. Seidel, 35 Y. Seiya, 39 
A. Semenov, 13 F. Sforza cc , 44 A. Sfyrla, 22 S.Z. Shalhout, 7 T. Shears, 27 P.F. Shepard, 45 M. Shimojima*, 53 
S. Shiraishi, 11 M. Shochet, 11 I. Shreyber, 34 A. Simonenko, 13 P. Sinervo, 31 A. Sissakian*, 13 K. Sliwa, 54 J.R. Smith, 7 
F.D. Snider, 15 A. Soha, 15 S. Somalwar, 50 V. Sorin, 4 P. Squillacioti, 15 M. Stancari, 15 M. Stanitzki, 59 
R. St. Denis, 19 B. Stelzer, 31 O. Stelzer-Chilton, 31 D. Stcntz, 36 J. Strologas, 35 G.L. Stryckcr, 32 Y. Sudo, 53 
A. Sukhanov, 16 I. Suslov, 13 K. Takemasa, 53 Y. Takeuchi, 53 J. Tang, 11 M. Tecchio, 32 P.K. Teng, 1 J. ThomV 5 
J. Thome, 10 GA. Thompson, 22 E. Thomson, 43 P. Ttito-Guzman, 29 S. Tkaczyk, 15 D. Toback, 51 S. Tokar, 12 
K. Tollcfson, 33 T. Tomura, 53 D. Tonclli, 15 S. Torre, 17 D. Torretta, 15 P. Totaro, 41 M. Trovato 66 , 44 Y. Tu, 43 
F. Ukegawa, 53 S. Uozumi, 25 A. Varganov, 32 F. Vazquez fe , 16 G. Velev, 15 C. Vellidis, 3 M. Vidal, 29 I. Vila, 9 
R. Vilar, 9 J. Vizan, 9 M. Vogel, 35 G. Volpi cc , 44 P. Wagner, 43 R.L. Wagner, 15 T. Wakisaka, 39 R. Wallny, 8 
S.M. Wang, 1 A. Warburton, 31 D. Waters, 28 M. Weinberger, 51 W.C Wester III, 15 B. Whitchouse, 54 D. Whiteson c , 43 
A.B. Wicklund, 2 E. Wickhmd, 15 S. Wilbur, 11 F. Wick, 24 H.H. Williams, 43 J.S. Wilson, 37 P. Wilson, 15 B.L. Winer, 37 
P. WittichV 5 S. Wolbers, 15 H. Wolfe, 37 T. Wright, 32 X. Wu, 18 Z. Wu, 5 K. Yamamoto, 39 J. Yamaoka, 14 

T. Yang, 15 U.K. Yang?, 11 Y.C. Yang, 25 W.-M. Yao, 26 G.P. Ych, 15 K. Yi" 1 , 15 J. Yoh, 15 K. Yorita, 56 
T. YoshidaJ', 39 G.B. Yu, 14 I. Yu, 25 S.S. Yu, 15 J.C. Yun, 15 A. Zanctti, 52 Y. Zcng, 14 and S. Zucchclli aa 61 

(CDF Collaboration^) 

1 Institute of Physics, Academia Sinica, Taipei, Taiwan 11529, Republic of China 
2 Argonne National Laboratory, Argonne, Illinois 60439, USA 
3 University of Athens, 157 71 Athens, Greece 
4 Institut de Fisica d'Altes Energies, ICREA, Universitat Autonoma de Barcelona, E-08193, Bellaterra (Barcelona), Spain 

5 Baylor University, Waco, Texas 76798, USA 
6 Istituto Nazionale di Fisica Nucleare Bologna, aa University of Bologna, 1-40127 Bologna, Italy 
7 University of California, Davis, Davis, California 95616, USA 
8 University of California, Los Angeles, Los Angeles, California 90024, USA 
Instituto de Fisica de Cantabria, CSIC-University of Cantabria, 39005 Santander, Spam 
10 Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA 
11 Enrico Fermi Institute, University of Chicago, Chicago, Illinois 60637, USA 
12 Comenius University, 842 48 Bratislava, Slovakia; Institute of Experimental Physics, 040 01 Kosice, Slovakia 
13 Joint Institute for Nuclear Research, RU-141980 Dubna, Russia 
14 Duke University, Durham, North Carolina 27708, USA 
15 Fermi National Accelerator Laboratory, Batavia, Illinois 60510, USA 
16 University of Florida, Gainesville, Florida 32611, USA 
17 Laboratori Nazionali di Frascati, Istituto Nazionale di Fisica Nucleare, 1-00044 Frascati, Italy 
18 University of Geneva, CH-1211 Geneva 4, Switzerland 
19 Glasgow University, Glasgow G12 8QQ, United Kingdom 
20 Harvard University, Cambridge, Massachusetts 02138, USA 
21 Division of High Energy Physics, Department of Physics, 
University of Helsinki and Helsinki Institute of Physics, FIN-00014, Helsinki, Finland 
22 University of Illinois, Urbana, Illinois 61801, USA 
28 The Johns Hopkins University, Baltimore, Maryland 21218, USA 
24 Institut fur Experimentelle Kernphysik, Karlsruhe Institute of Technology, D-76131 Karlsruhe, Germany 
25 Center for High Energy Physics: Kyungpook National University, 
Daegu 702-701, Korea; Seoul National University, Seoul 151-742, 
Korea; Sungkyunkwan University, Suwon 440-746, 
Korea; Korea Institute of Science and Technology Information, 
Daejeon 305-806, Korea; Chonnam National University, Gwangju 500-757, 
Korea; Chonbuk National University, Jeonju 561-756, Korea 
26 Ernest Orlando Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA 
27 University of Liverpool, Liverpool L69 7ZE, United Kingdom 
28 University College London, London WC1E 6BT, United Kingdom 
29 Centro de Investigaciones Energeticas Medioambientales y Tecnologicas, E-28040 Madrid, Spain 
30 Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA 
31 Institute of Particle Physics: McGill University, Montreal, Quebec, 
Canada H3A 2T8; Simon Fraser University, Burnaby, British Columbia, 
Canada V5A 1S6; University of Toronto, Toronto, Ontario, 
Canada M5S 1A7; and TRIUMF, Vancouver, British Columbia, Canada V6T 2A3 
32 University of Michigan, Ann Arbor, Michigan 48109, USA 
33 Michigan State University, East Lansing, Michigan 48824, USA 
34 Institution for Theoretical and Experimental Physics, ITEP, Moscow 117259, Russia 



3 



35 University of New Mexico, Albuquerque, New Mexico 87131, USA 
36 Northwestern University, Evanston, Illinois 60208, USA 
S7 The Ohio State University, Columbus, Ohio 43210, USA 
38 Okayama University, Okayama 700-8530, Japan 
Osaka City University, Osaka 588, Japan 
40 University of Oxford, Oxford 0X1 3RH, United Kingdom 
41 Istituto Nazionale di Fisica Nucleare, Sezione di Padova-Trento, bb University of Padova, 1-35131 Padova, Italy 
4 ' 2 LPNHE, Unwersite Pierre et Marie Curie/IN2P3-CNRS, UMR7585, Paris, F-75252 France 
43 University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA 

44 Istituto Nazionale di Fisica Nucleare Pisa, cc University of Pisa, 

dd University of Siena and ee Scuola Normale Superiore, 1-56127 Pisa, Italy 

45 University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA 

46 Purdue University, West Lafayette, Indiana 47907, USA 
47 University of Rochester, Rochester, New York 14627, USA 
48 The Rockefeller University, New York, New York 10065, USA 
49 Istituto Nazionale di Fisica Nucleare, Sezione di Roma 1, 

Sapienza Universitd di Roma, 1-00185 Roma, Italy 
50 Rutgers University, Piscataway, New Jersey 08855, USA 
51 Texas A&M University, College Station, Texas 77843, USA 

52 Istituto Nazionale di Fisica Nucleare Trieste/ U dine, 
1-34100 Trieste, 33 University of Udine, 1-33100 Udme, Italy 
53 University of Tsukuba, Tsukuba, Ibaraki 305, Japan 
54 Tufts University, Medford, Massachusetts 02155, USA 
55 University of Virginia, Charlottesville, Virginia 22906, USA 
5b Waseda University, Tokyo 169, Japan 
57 Wayne State University, Detroit, Michigan 48201, USA 
University of Wisconsin, Madison, Wisconsin 53706, USA 

59 Yale University, New Haven, Connecticut 06520, USA 

60 University of Virginia, Charlottesville, VA 22906, USA 

61 Istituto Nazionale di Fisica Nucleare Bologna, aa University of Bologna, 1-40127 Bologna, Italy 

(Dated: September 13, 2011) 

We present a study of B® decays to the CP-odd final state J/ipfo (980) with J/0 — » fi + fi~ and 
/o(980) — > 7r + 7r~. Using pp collision data with an integrated luminosity of 3.8 fb collected by the 
CDF II detector at the Tevatron we measure a B° lifetime of t(B° -> J/0/o(98O)) = l.TOj^ii (stat)± 
0.03(syst) ps. This is the first measurement of the B® lifetime in a decay to a CP eigenstate and 
corresponds in the standard model to the lifetime of the heavy B® eigenstate. We also measure the 
product of branching fractions of B® J/0 /o (980) and /o(980) — > tt + it~ relative to the product of 
branching fractions of B° -> J/00 and -» K + K~ to be R fo/4> = 0.257 ± 0.020(stat) ± 0.014(syst), 
which is the most precise determination of this quantity to date. 
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In the standard model, the mass and flavor eigenstates 
of the B® meson are not identical. This gives rise to 
particle - anti-particle oscillations which proceed in 
the standard model through second order weak interac- 
tion processes, and whose phenomenology depends on 
the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing 
matrix. The time (t) evolution of B® mesons is approxi- 
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mately governed by the Schrodinger equation 



dt 



\B° s (t)) 
\B M 



M s 



\B° s (t)) 
\B°M 



(1) 



where M s and t s are mass and decay rate symmetric 
2x2 matrices. Diagonalization of M s 
mass eigenstates 



\B" S 
\B[ 



Hi 



= P \B l s 
= P \B C S 



l\B° s ) 

q \B°s) 



\T S leads to 

(2) 
(3) 



with distinct masses (M^, ) and distinct decay rates 
(T^, r, ), where p and q are complex numbers satisfying 



\P\ 2 



1. An important feature of the B Q S system is 



the non-zero matrix element Y\ 2 representing the partial 
width of and B° s decays to common final states which 
translates into a non-zero decay width difference Al? s of 
the two mass eigenstates through the relation 



r 



H 



(4) 



where <fi s = arg(— APf 2 /rf 2 ). The phase <f> 8 describes CP 
violation in B® mixing. In the standard model <p s is pre- 
dicted to be 0.22° ± 0.06° @, Q. The small value of the 
phase (j) s causes the mass and CP eigenstates to coin- 
cide to a good approximation. Thus the measurement 
of the lifetime in a CP eigenstate provides directly the 
lifetime of the corresponding mass eigenstate. If new 
physics is present, it could enhance 4> s to large values, a 
scenario which is not excluded by current experimental 
constraints. In such a case the correspondence between 
mass and CP eigenstates does not hold anymore and the 
measured lifetime will correspond to the weighted av- 
erage of the lifetimes of the two mass eigenstates with 
weights dependent on the size of the CP violating phase 
4>s 0] ■ Thus a measurement of the B® lifetime in a final 
state which is a CP eigenstate provides, in combination 
with other measurements, valuable information on the 
decay width difference Ar s and the CP violation in P° 
mixing. 

One of the most powerful measurements to constrain 
a new physics contribution to the phase <p s is the mea- 
surement of CP violation in the decay B® — > J/tp4> with 
<f> — > K + K~ . The decay P° — > J/ijjcj) has a mixture of the 
CP-even and -odd components in the final state and an 
angular analysis is needed to separate them In the 
standard model, CP violation in the decay B® — > J/tpcf) 
is given by /3 S = &Yg[(-V ts V* b ) / {V cs V* b )]. New physics 
effects in P" mixing would shift (f> s and — 2/3 s from the 
standard model value by the same amount. A sufficiently 
copious B® —> J/ipfo signal with /q — > ir + ir~, where /o 
stands for Jo (980), and P° flavor identified at production 
can be used to measure /3 S without the need of an angu- 
lar analysis Q as J/ipfo is a pure CP-odd final state. 
Since the P° is a spin particle and the decay prod- 
ucts J/tp and /o have quantum numbers J PC = 1 and 
++ , respectively, the final state has an orbital angu- 
lar momentum of L = 1 leading to a CP eigenvalue of 



(— 1) L = — 1. Further interest in the decay P° — > J/ipfo 
arises from its possible contribution to an 5-wave com- 
ponent in the B° s — > J/ipK + K~ decay if the /o decays 
to K + K~ . This contribution could help to resolve an 
ambiguity in the Ar s and f3 s values determined in the 
P" — > J/4>4> analyses. Because it was neglected in the first 
tagged B® — > J/ipcj) results @, H| , each of which showed an 
approximately 1.5 a deviation from the standard model, 
it was argued that the omission may significantly bias the 
results [UH3]. However, using the formalism in Ref. [Tl| . 
the latest preliminary CDF measurement [l2[ has shown 
that the S-w&ve interference effect is negligible at the 
current level of precision. 

In Refs. pi, 3] the decay width difference in the stan- 
dard model is predicted to be Arf M = (0.087 ± 0.021) 
ps _1 and the ratio of the average P° lifetime, t s = 
2/(rf +rf ), to the B° lifetime, r d , to be 0.996 < T s /r d < 
1. Using these predictions in the relations 
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7h 
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2' 
1 



-Ar, 



-Ar e 



(5) 
(G) 



where T s = 1/t s , together with the world average B° 
lifetime, r d = (1.525 ± 0.009) ps Q, we find the 
theoretically-derived values = (1.630 ± 0.030) ps and 
7f = (1.427 ±0.023) ps. 

While no direct measurements of Pi? lifetimes in decays 
to pure CP eigenstates are available, various experimen- 
tal results allow for the determination of the lifetimes 
of the two mass eigenstates. Measurements sensitive to 
these lifetimes are the angular analysis of B° s — > J/ip(f> 
decays and the branching fraction of B® — > Pi*' + D B , 
which can be complemented by measurements of the P° 
lifetime in flavor specific final states. The combination 
of available measurements yields t^ 1 — (1.544± 0.041) ps 
and t s l = (1.4071^026) PS [l!- From CDF measurements 
we can infer the two lifetimes from the result of the angu- 
lar analysis of P° — > J/tp(f> decays. The latest preliminary 
result [12| , that is not yet included in the above average, 
yields rf = (1.622 ± 0.068) ps and t s l = (1.446 ± 0.035) 
ps assuming standard model CP violation. 

Compared to measurements using P" — > J/?p(j) de- 
cays, lifetime and future CP violation measurements in 
the P" — > J/ipfo decay suffer from a lower branching 
fraction. Based on a comparison to Df meson decays 
Ref. [9J] makes a prediction for the branching fraction of 
Pi? — > J/0/o decay relative to the P° — > J/ipcf> decay, 



P 



fo/<P 



B(B° S -> J/0/ o ) B(f 



B{B° S Jftpcp) B{<p -> K+K- 



(7) 



to be approximately 0.2. The CLEO experiment es- 
timates Rfg/4, — 0.42 ± 0.11 from a measurement of 
scmilcptonic decays [l6|. A theoretical prediction 
based on QCD factorization yields a range of P/ /«/, be- 
tween 0.2 and 0.5 [l7|. With the world average branching 
fraction for the P° ->■ J/ip(f> decay of (1.3 ± 0.4) x 10" 3 
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and the branching fraction of fo —> tt + it in the re gion 
between 0.5-0.8, predictions of B(B° S -> J/0/ o ) [Hill 
translate into a wide range of i?/ o /0 values of approxi- 
mately 0.1-0.5. 

The first experimental search was performed by the 
Belle experiment (20j . Their preliminary result did not 
yield a signal and they extract an upper limit on the 
branching fraction of B(B° -> J/ipf )B(fo -> 7r+7r~) < 
1.63 x 10~ 4 at90%C.L. Recently the LHCb experi- 
ment reported the first observation of the decay B® — > 
J/ipfo |2l| with a relative branching fraction of Rf /$ — 
0.252j£^(stat) to.°33(syst). Shortly after the LHCb re- 
sult was presented, the Belle collaboration announced 
their result of an updated analysis using 121.4 fb _1 of 
T(55) data [22| . They observe a significant B° s — > J/ipfa 
signal and measure B(B° — > J/ipfo)B(fa — > tt + tt~) = 

(1.16±g:?| ±8^17 iolfi) x 10_4 > where the first uncertainty 
is statistical, the second systematic, and the third one 

is an uncertainty on the number of produced B s B s 
pairs. Using their preliminary measurement of the B® — > 
J/ip(j) branching fraction [23} , and assuming that the un- 
certainty on the number of produced B s ' B s pairs 
is fully correlated for the two measurements, this trans- 
lates into R fa/(j> = 0.206±g;g!i(stat) ± 0.052(syst). A 
preliminary measurement of the DO experiment yields 
R fo/<t> = 0.210 ± 0.032(stat) ± 0.036(syst) [p. 

In this paper we present a measurement of the ratio 
Rf /4> of the branching fraction of the B® — > J/ipfo decay 
relative to the J5° — > J/iptp decay and the first measure- 
ment of the Bg lifetime in a decay to a pure CP eigen- 
state. We use data collected by the CDF II detector from 
February 2002 until October 2008. The data correspond 
to an integrated luminosity of 3.8 fb . 

This paper is organized as follows: In Sec. |TT] we de- 
scribe the CDF II detector together with the online data 
selection, followed by the candidate selection in Sec. Mil 
Section IIVI describes details of the measurement of the 
ratio Rf /4> of branching fractions of the B° s — > J/tpfo 
decay relative to the B® —> J/ip(f> decay while Sec. |V] dis- 
cusses the lifetime measurement. We finish with a short 
discussion of the results and conclusions in Sec. IVIl 



II. CDF II DETECTOR AND TRIGGER 

Among the components of the CDF II detector [25| 
the tracking and muon detection systems are most rel- 
evant for this analysis. The tracking system lies within 
a uniform, axial magnetic field of 1.4 T strength. The 
inner tracking volume hosts 7 layers of double-sided sil- 
icon micro-strip detectors up to a radius of 28 cm [26| . 
An additional layer of single-sided silicon is mounted di- 
rectly on the beam-pipe at a radius of 1.5 cm, provid- 
ing an excellent resolution of the impact parameter do, 
defined as the distance of closest approach of the track 
to the interaction point in the transverse plane. The 
silicon tracker provides a pseudorapidity coverage up to 



\r)\ < 2.0. The remainder of the tracking volume up to a 
radius of 137 cm is occupied by an open-cell drift cham- 
ber (27|- The drift chamber provides up to 96 measure- 
ments along the track with half of them being axial and 
other half stereo. Tracks with |^| < 1.0 pass the full 
radial extent of the drift chamber. The integrated track- 
ing system achieves a transverse momentum resolution 
of (j{pt)/Pt ~ 0.07% (GeV/c) -1 and an impact param- 
eter resolution of a (do) « 35 /xm for tracks with a trans- 
verse momentum greater than 2 GeV/c. The tracking 
system is surrounded by electromagnetic and hadronic 
calorimeters, which cover the full pseudorapidity range 
of the tracking system [28M3l1 ] . We detect muons in three 
sets of multi-wire drift chambers. The central muon de- 
tector has a pseudorapidity coverage of \rj\ < 0.6 [32[ 
and the calorimeters in front of it provide about 5.5 in- 
teraction lengths of material. The minimum transverse 
momentum to reach this set of muon chambers is about 
1.4 GeV/c. The second set of chambers covers the same 
range in rj, but is located behind an additional 60 cm of 
steel absorber, which corresponds to about 3 interaction 
lengths. It has a higher transverse momentum threshold 
of 2 GeV/c, but provides a cleaner muon identification. 
The third set of muon detectors extends the coverage to 
a region of 0.6 < \r)\ < 1.0 and is shielded by about 6 
interaction lengths of material. 

A three-level trigger system is used for the online event 
selection. The trigger component most important for this 
analysis is the extremely fast tracker (XFT) (33[, which 
at the first level groups hits from the drift chamber into 
tracks in the plane transverse to the beamline. Candidate 
events containing J/tp decays are selected by a 

dimuon trigger [25J which requires two tracks of opposite 
charge found by the XFT that match to track segments 
in the muon chambers and have a dimuon invariant mass 
in the range 2.7 to 4.0 GeV/c 2 . 



III. RECONSTRUCTION AND CANDIDATE 
SELECTION 

A. Reconstruction 

In the offline reconstruction we first combine two muon 
candidates of opposite charge to form a J/ip candidate. 
We consider all tracks that can be matched to a track 
segment in the muon detectors as muon candidates. The 
J/ip candidate is subject to a kinematic fit with a vertex 
constraint. We then combine the J/ip candidate with two 
other oppositely charged tracks that are assumed to be 
pions and have an invariant mass between 0.85 and 1.2 
GeV/c 2 to form a 6 S ° -> J/ipfo candidate. In the final 
step a kinematic fit of the B® — > J/ipfo candidate is per- 
formed. In this fit we constrain all four tracks to originate 
from a common vertex, and the two muons forming the 
J/0 are constrained to have an invariant mass equal to 
the world average J/0 mass [b"i| . In a similar way we also 
reconstruct B° s — > J/ip<f) candidates using pairs of tracks 
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of opposite charge assumed to be kaons and having an 
invariant mass between 1.009 and 1.029 GeV/c 2 . Dur- 
ing the reconstruction we place minimal requirements on 
the track quality, the quality of the kinematic fit, and 
the transverse momentum of the B® candidate to ensure 
high quality measurements of properties for each candi- 
date. For the branching fraction measurement we add 
a requirement which aims at removing a large fraction 
of short-lived background. We require the decay time of 
the B® candidate in its own rest frame, the proper decay 
time, to be larger than three times its uncertainty. This 
criterion is not imposed in the lifetime analysis since it 
would bias the lifetime distribution. The proper decay 
time is determined by the expression 

C-Pt 

where L xy is the flight distance projected onto the B° s 
momentum in the plane transverse to the beamline, Pt 
is the transverse momentum of the given candidate, and 
m(B a ) is the reconstructed mass of the B° candidate. 
The uncertainty on the proper decay time t is estimated 
for each candidate by propagating track parameter and 
primary vertex uncertainties into an uncertainty on L xy . 
The proper decay time resolution is typically of the order 
ofO.lps. 



B. Selection 

The selection is performed usin g a neural network 
based on the NEUROBAYES package [34], Hf| • The neural 
network combines several input variables to form a single 
output variable on which the selection is performed. The 
transformation from the multidimensional space of input 
variables to the single output variable is chosen during 
a training phase such that it maximizes the separation 
between signal and background distributions. For each 
of the two measurements presented in this paper we use 
a specialized neural network. For the training we need 
two sets of events with a known classification of signal 
or background. For the signal sample we use simulated 
events. We generate the kinematic distributions of B° 
mesons according to the measured 6-hadron momentum 
distribution. The decay of the generated B® particles 
into the Jfyfo final state is simulated using the evtgen 
package [36| . Each event is passed through the standard 
CDF II detector simulation, based on the GEANT3 pack- 
age [13, Hll . The simulated events are reconstructed with 
the same reconstruction software as real data events. The 
background sample is taken from data using candidates 
with the J/ipir + ir~ invariant mass above the B® signal 
peak, where only combinatorial background events con- 
tribute. Because the requirement on the proper decay 
time significance efficiently suppresses background events 
in the branching ratio measurement, we use an enlarged 
sideband region of 5.45 to 5.55 GeV/c 2 in this analysis, 



compared to an invariant mass range from 5.45 to 5.475 
GeV/c 2 for the lifetime measurement. 

For the branching fraction measurement, the inputs to 
the neural network, ordered by the importance of their 
contribution to the discrimination power, are the trans- 
verse momentum of the /o, the \ 2 of the kinematic fit of 
the B® candidate using information in the plane trans- 
verse to the beamline, the proper decay time of the B® 
candidate, the quality of the kinematic fit of the B® can- 
didate, the helicity angle of the positive pion, the trans- 
verse momentum of the B® candidate, the quality of the 
kinematic fit of the two pions with a common vertex con- 
straint, the helicity angle of the positive muon, and the 
quality of the kinematic fit of the two muons with com- 
mon vertex constraint. The helicity angle of the muon 
(pion) is defined as the angle between the three momenta 
of the muon (pion) and B° candidate measured in the rest 
frame of the J/ip (fo). For the selection of B® — > J/ip(f> de- 
cays we use the same neural network without retraining 
and simply replace /o variables by <fi variables and pions 
by kaons. 

For the lifetime measurement we modify the list of in- 
puts by removing the proper decay time. We also do not 
use the helicity angles as they provide almost no addi- 
tional separation power on the selected sample. Since we 
are not concerned about a precise efficiency determina- 
tion for the lifetime measurement, we add the following 
inputs: the invariant mass of the two pions, the likelihood 
based identification information for muons [39| . and the 
invariant mass of the muon pair. The muon identifica- 
tion is based on the matching of tracks from the tracking 
system to track segments in the muon system, energy 
deposition in the electromagnetic and hadronic calorime- 
ters, and isolation of the track. The isolation is defined as 
the transverse momentum carried by the muon candidate 
over the scalar sum of transverse momenta of all tracks in 
a cone of AR = ^{A<t>) 2 + (Ar/) 2 < 0.4, where Acj) (At)) 
is the difference in azimuthal angle (pseudorapidity) of 
the muon candidate and the track. There is no signif- 
icant change in the importance ordering of the inputs. 
The invariant mass of the pion pair becomes the second 
most important input, the likelihood based identification 
of the two muon candidates is ranked fourth and sixth in 
the importance list, and the muon pair invariant mass is 
the least important input. 

For the branching fraction measurement we select the 
threshold on the neural network output by maximizing 
e/(2.5 + ^/Nb) .40], where e is the reconstruction effi- 
ciency for B® — > J/ipfo decays and Nb is the number 
of background events estimated from the J/07r + 7r~ mass 
sideband. The invariant mass distributions of selected 
B® — > J/ipfo and B® — > J/tp(f> candidates are shown in 
Figs. [1] and O A clear signal at around 5.36 GeV/c 2 is 
visible in both mass distributions. 

For the lifetime measurement we use simulated exper- 
iments to determine the optimal neural network output 
requirement. We select a value that minimizes the sta- 
tistical uncertainty of the measured lifetime. We scan 



7 




CM 

o 


0.045 r 


> 


- 


CD 


0.040^ 


m 
i 


0.035 L 


CD 




Q_ 


0.030 1 






'3 


0.025 1 


CO 




_Q 

O 


0.020 1 


i 

Q. 


0.015I 




0.010 




0.005 1 




0.000^ 




^ B"^ JA|/ K 

Other B° 

1 bS-^j/v* 

B° s -> JA|< K + K 
|=| Other B° 



5.4 



5.5 



M(J/ V 7i + 7t ) [GeV/c' 



M(JA|m + 7f) [GeV/c 2 



FIG. 1: The invariant mass distribution of — > J/ipfo can- 
didates selected for the branching fraction measurement. 
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FIG. 2: The invariant mass distribution of — > J/ip(l> candi- 
dates selected for the branching fraction measurement. 



a wide range of neural network output values and for 
each requirement we simulate an ensemble of experiments 
with a lifetime of 1.63 ps, where the number of signal 
and background events as well as the background distri- 
butions are simulated according to data. For a broad 
range of selection requirements we observe the same un- 
certainty within a few percent. Our final requirement on 
the network output is chosen from the central region of 
this broad range of equivalent options. 



FIG. 3: (color online) Stacked histogram of physics back- 
grounds to B° —> Jfipfo derived from simulation using the 
selection for the branching fraction measurement. The verti- 
cal line indicates the location of the world average B° mass. 



we use the selection from the branching fraction measure- 
ment. While several physics backgrounds appear in the 
J/i/jtt + tt~ mass spectrum, none contributes significantly 
under the B° s peak. The most prominent physics back- 
grounds are B° — > J/ipK*° with K*° -> K + ir~ , where 
A"* stands for K* (892)°, and B° ->■ J/V-tt+tt" . In the 
first case the kaon is mis-reconstructed as a pion and 
gives rise to a large fraction of the structure seen be- 
low 5.22 GeV/c 2 , while the second one is correctly re- 
constructed and produces the narrow peak at approxi- 
mately 5.28 GeV/c 2 . Another possible physics type of 
background would consist of properly reconstructed B + 
combined with a random track. This type of background 
would contribute only to higher masses with a thresh- 
old above the signal. As we do not find evidence 
of such background in Ref. [l3[ which is more sensitive 
we conclude that this kind of background is also negligi- 
ble here. The stacked histogram of physics backgrounds 
derived from simulation is shown in Fig. [3] From this 
study we conclude that the main physics background 
that has to be included as a separate component in a 
fit to the mass spectrum above 5.26 GeV/c 2 stems from 
decays of B° —> J/-i/j7t + it~ . It is properly reconstructed 
and therefore simple to parametrize. All other physics 
backgrounds are negligible. 



IV. BRANCHING FRACTION MEASUREMENT 



C. Physics backgrounds 

We study possible physics backgrounds using simu- 
lated events with all 6-hadrons produced and decayed 
inclusively to final states containing a J/ip. For this study 



In this Section we describe details of the branching 
fraction measurement. These involve the maximum like- 
lihood fit to extract the number of signal events, the effi- 
ciency estimation, and the systematic uncertainties. We 
conclude this Section with the result for the ratio -R/ /</> of 
branching fractions between B® — > J/ipfo and B® J/ip(/> 
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decays. 



A. Fit description 

We use an unbinned extended maximum likelihood fit 
of the invariant mass to extract the number of B® decays 
in our samples. In order to avoid the need for modeling 
most of the physics background, we restrict the fit to 
the mass range from 5.26 GeV/c 2 to 5.5 GeV/c 2 . The 
likelihood is 

N 

£ = I1 [N s -P s (mi) + N cb -P cb (mi)+ 
i=i 

fpb ■ N s ■ P p b(mi) + N B o ■ P B o(mi)] ■ 

e -(N s +N ob +N s -f pb +N B o) ( fg\ 

where m, is the invariant mass of the i-th candidate and 
N is the total number of candidates in the sample. The 
fit components are denoted by the subscripts s for sig- 
nal, cb for combinatorial background, pb for physics back- 
ground, and B° for B° — > J/i(jtt + tt~ background. The 
yields of the components are given by N s , N c b, N s ■ f p b, 
and N B o , and their probability density functions (PDFs) 
by P s ,cb, P b,B°{' m i)j respectively. The physics background 
yield is parametrized relative to the signal yield via the 
ratio fpb to allow constraining it by other measurements 
in the B° s -> J/if)<f> fit. 

The signal PDF P s (mj) is parametrized by a sum of 
two Gaussian functions with a common mean. The rela- 
tive size of the two Gaussians and their widths are deter- 
mined from simulated events. Approximately 82% of the 
Bg —> J/ipfo decays are contained in a narrower Gaus- 
sian with width of 9.4 MeV/c 2 . The broader Gaussian 
has width of 18.4 MeV/c 2 . In the case of Bg — > J/ipcf), the 
narrow Gaussian with a width of 7.2 MeV/c 2 accounts 
for 79% of the signal, with the rest of the events having 
a width of 13.3 MeV/c 2 . To take into account possible 
differences between simulation and data, we multiply all 
widths by a scaling parameter S m . Because of kinematic 
differences between f —> ir + n~ and <f) — > K + K~ we use 
independent scale factors for both modes. In the fits all 
parameters of the PDF are fixed except for the scaling 
parameter S m . In addition the mean of the Gaussians is 
allowed to float in the J/ipK + K~ fit. Doing so we obtain 
a value that is consistent with the world average B® mass 
[l4j . For the J/ipir + iT~ fit we fix the position of the sig- 
nal to the value determined in the fit to the J/ipK + K~ 
candidates. 

The combinatorial background PDF P c b(mi) is 
parametrized using a linear function. In both fits we 
leave its slope floating. In each of the two fits there is one 
physics background. In the case of the J/07r + 7r~ spec- 
trum, the physics background describes properly recon- 
structed B° —> J/tpTT + TT~ decays using a shape identical 
to the B® sig nal and position fixed to the world average 
B° mass [14|. The number of B° events Ngo is left free 



in the fit. For the J/il)K + K~ fit, we have a contribution 
from B° — > J/ipK*° decays where the pion from the A"* 
decay is mis-reconstructed as a kaon. This contribution 
peaks at a mass of approximately 5.36 GeV/c 2 with an 
asymmetric tail towards larger masses. The shape itself 
is parametrized by a sum of a Gaussian function and an 
exponential function convolved with a Gaussian. The 
parameters are derived from simulated B° — > J/ipK*° 
events. The normalization of this component relative to 
the signal is fixed to f pb = (3.04 ± 0.99) x 10~ 2 , which is 
derived from the CDF Run I measurement of the ratio of 
cross section times branching fraction for B® — > J/0</> and 
B° — > J/ijj K*° decays I4ll. the world average branching 
fractions for <j) and K*° [14| , and the ratio of reconstruc- 
tion efficiencies obtained from simulation. 

The fit determines a yield of 502 ± 37 B° s ->• J/0/ o 
events and 2302 ± 49 B® —> J/tp4> events, where the un- 
certainties are statistical only. The number of B° back- 
ground events in the J/07r + 7r~ fit is 160 ± 30. 



B. Efficiency 

To extract the ratio of branching fractions we need to 
estimate the relative efficiency for reconstruction of B° — > 
Jfipfo with f -> tt+tt" and B° -s- J/ipcj) with -S- A+A~ 
decays, e rei = e(B° s -> J/V>0)/e(B° J/i/jf ). We esti- 
mate the efficiency using simulated events in which we 
generate a single B® meson per event. The B° a meson 
then decays with equal probabilities to B® — > J/V'/o or 
B® — > J/ip4> final states with exclusive J/ip — > 
cf> — > A + A~, and /o — > 7r + 7r~. Generated events are 
then processed through a detailed detector simulation 
and the offline reconstruction software used to recon- 
struct data. In both cases angular and decay time dis- 
tributions are generated assuming no CP violation and 
parameters taken from the preliminary result of the an- 
gular distributions analysis 12]: r = 1.529 ± 0.028 ps, 
Ar = 0.075 ± 0.036 ps" 1 , L4 | 2 = 0.524 ± 0.020, and 
|^4||| 2 = 0.231 ± 0.021. As a strong phase between A 
and A|| is not measured we use the world average value 
from B° -> J/tpK* decays of_0 = -2.86 ± 0.11 [3 as 
a reasonable approximation [42| . An additional peculiar- 
ity of the B® — > Jfyfo decay is the unusual mass shape 
of the ,fo meson. It is modeled using a Flatte distribu- 
tion [431 ] with input parameters measured by the BES 
experiment [H| to be m = 965 ± 8 ± 6 MeV/c 2 , g n = 
165 ± 10 ± 15 MeV/c 2 , and g K /g v = 4.21 ± 0.25 ± 0.21, 
where the errors are statistical and systematic, respec- 
tively. The <j) meson mass distribution is modeled using 
a relativistic Breit-Wigner distribution with world aver- 
age values for its parameters [3]. With these inputs to 
the simulation we find e re i = 1.178, which accounts for 
the (f> an d fo mass window selection requirements. 
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C. Systematic uncertainties 

We investigate several sources of systematic uncertain- 
ties. They can be broadly separated into two classes: one 
dealing with assumptions made in the fits that may af- 
fect yields, and the other related to assumptions in the 
efficiency estimation. In the first class we estimate uncer- 
tainties by refitting data with a modified assumption and 
taking the difference with respect to the original value 
as an uncertainty. For the second class we recalculate 
the efficiency with a modified assumption and take the 
difference with respect to the default efficiency as an un- 
certainty unless specified otherwise. The summary of as- 
signed uncertainties is given in Table HI 

For the yield of B® — > J/ip'P we investigate the ef- 
fect of the assumption on the combinatorial background 
shape, the limited knowledge of mis-reconstructed B° — > 
J/ipK*° decays and the shape of the signal PDF. The un- 
certainty due to the shape of combinatorial background is 
estimated by changing from the first order polynomial to 
a constant or a second order polynomial. For the physics 
background we vary the normalization of the component 
in the fit and use an alternative shape determined by 
varying the momentum distribution and the decay am- 
plitudes of B° — > J/ipK*° in simulation. Finally, to esti- 
mate the effect of the signal PDF parametrization we use 
an alternative model with a single Gaussian rather than 
two Gaussian functions and an alternative shape from 
simulation, where we vary the momentum distribution of 
the produced B® mesons and the decay amplitudes of the 
B° s -> J/00 decay. 

To estimate the uncertainty on the B° — ► J/0/o yield 
we follow a procedure similar to that for B® — ¥ J/ip<fi an d 
conservatively treat the systematic effects as independent 
between the two modes in the calculation of Rf / t f > . For 
the sensitivity to the parametrization of the combinato- 
rial background we switch to a second order polynomial 
or a constant as alternative parametrization. For the 
shape of the signal PDF we use two alternatives, one 
with a single Gaussian function instead of two and an- 
other one with two Gaussians, but varying the momen- 
tum distribution in simulation. We also vary the position 
of the B® signal within the uncertainty determined in the 
J/ipK + K- fit. 

The systematic uncertainty on the relative efficiency 
stems from the statistics of simulation, an imperfect 
knowledge of the momentum distribution, physics pa- 
rameters of decays like lifetimes or decay amplitudes, 
and differences in the efficiencies of the online selection 
of events. To estimate the effect of the imperfect knowl- 
edge of the momentum distribution we vary the momen- 
tum distribution of B® mesons in the simulation. The 
physics parameters entering the simulation are grouped 
into three categories, those defining the /o mass shape, 
the ones determining decay amplitudes in B® — > J/ip(f> de- 
cays, and those affecting the lifetimes of the two B® mass 
eigenstates. In the first two cases we vary each parame- 
ter independently and add all changes in the efficiency in 



TABLE I: The summary of assigned systematic uncertainties 
for the branching fraction measurement. 



Source 


J/0</> yield J/0/o yield 




Combinatorial bckg. 


34 


16 


— 


Physics bckg. 


13 






iviabo lehoiuLiuii 


32 


7.9 




_B° mass 




0.1 




Total 


49 


18 




MC statistics 






0.012 


Momentum distribution 






0.011 


Decay parameters 






0.033 


Trigger composition 






0.016 


Total 






0.040 



quadrature. For the last case we vary the mean lifetime 
t and the decay width difference Ar simultaneously and 
take the largest variation as the uncertainty. We add the 
uncertainty from the third class in quadrature with all 
others to obtain the uncertainty due to the parameters 
describing the particle decays. The last effect deals with 
how events are selected during data taking. The CDF 
trigger has several different sets of requirements for the 
selection of events. The ones used in this analysis can be 
broadly sorted into three classes depending on momen- 
tum thresholds and which subdetectors detected muons. 
The fraction of events for each different class varies de- 
pending on the instantaneous luminosity, which is not 
simulated. To estimate the size of a possible effect we 
calculate the efficiency for each class separately and take 
half of the largest difference as the uncertainty. 

To obtain the total uncertainty we add all partial un- 
certainties in quadrature. In total we assigned 49 events 
(2.1%) as the systematic uncertainty on the B® — > J/00 
yield, 18 events (3.6%) on the B° — > J/tjjfo yield, and 
0.040 (3.4%) on the relative efficiency e re ;. A summary 
of the systematic uncertainties in the branching ratio is 
provided in Table HI 



D. Branching fraction result 

From the fit we find 502 ± 37(stat) ± 18(syst) B° -> 
J/0/o signal events and 2302 ± 49(stat) ± 49(syst) 5° 
J/00 events. The projections of the fits for B® — > J/4>fo 
and B® — > Jj^4> are shown in Fig. [4] and Fig. [5l respec- 
tively. 

In order to check our interpretation of the signal in 
the J/ipir + Tr~ distribution being due to the B° -> J/0/ o 
decays we show the invariant mass distribution of the 
pions for B° s signal data in Fig. |6l To obtain the dis- 
tribution of B® signal we fit the J/07r + 7r~ mass distri- 
bution in the range 5.26 to 5.45 GeV/c 2 for each bin in 
7r + 7r~ mass and report the B Q S signal yield as a function of 
7r + 7r~ mass. We fit the dipion mass distribution using the 
Flatte parametrization. The fit probability is 23.4% and 
the obtained parameters, too — 989.6±9.9(stat) MeV/c 2 , 

= 141±19(stat) MeV/c 2 , and g K /g n = 2.3±1.3(stat), 
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FIG. 4: (color online) Projection of the fit of the B° — > J/ipfo 
decay mode. The dashed line (blue) shows the contribution 
from combinatorial background. 
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FIG. 6: The dipion invariant mass distribution after sideband 
subtraction with fit projection overlaid. The fit uses a Flatte 
distribution with all parameters floating. 
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FIG. 5: (color online) Projection of the fit of the B® — > J/ijxj) 
decay mode. The dashed line (blue) and filled area (green) 
show the contributions from combinatorial background and 
B° — > J/ipK*°, respectively. 



are in reasonable agreement with the ones measured by 
the BES collaboration [Hj]. In Figs. [7] and |8] we show 
the positive muon and pion helicity angle distributions, 
obtained in an analogous way to the invariant mass dis- 
tribution of pion pairs. Those are corrected for relative 
efficiencies in the different helicity bins and compared to 
the theoretical expectation for aB s °-> J/ipfo signal. We 
use a x 2 test to evaluate the agreement between data and 
theoretical expectation. For the distribution of cos^^) 
we obtain x 2 / n df = 7.9/20, which corresponds to 99% 
probability. Similarly for cos^^) the x 2 / n df is 15/20, 




cos(e r ) 



FIG. 7: Normalized helicity angle distribution for the positive 
muon corrected for relative efficiency. The line shows the 
expectation for a B® — > J/ipfo decay. 



giving 78% probability. Since the dipion mass as well 
as the angular distributions are consistent with expecta- 
tions, we interpret our signal as coming solely from the 
£>!? — > J/ipfo decays. On the other hand, as we use a 
dipion mass window from 0.85 to 1.2 GeV/c 2 , we cannot 
exclude contributions from other higher mass states to 
our signal with present statistics. 

Finally, we obtain the ratio of branching fractions 



R 



fo/4> 



B(B° S -> J/#) B{4> -> K+K-) ' 
0.257 ± 0.020(stat) ± 0.014(syst), 



(10) 



where corrections for events with an f or ef> mass outside 
the ranges selected in this analysis are taken into account. 
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FIG. 8: Normalized helicity angle distribution for the posi- 
tive pion corrected for relative efficiency. The line shows the 
expectation for a — > Jfipfo decay. 



V. LIFETIME MEASUREMENT 

In this Section we discuss the details of the lifetime 
measurement. We describe the maximum likelihood fit, 
estimate the systematic uncertainties, and present the 
result of the lifetime measurement. 



A. Fit description 

To extract the B® lifetime we use a maximum likeli- 
hood fit. The fit uses three variables: the invariant mass 
m», the decay time U, and the decay time uncertainty au 
of each candidate. To exclude B° — > J/ipir + n~ decays we 
use only candidates with an invariant mass greater than 
5.3 GcV/c 2 in the fit. 

The components in the fit are B® signal and combina- 
torial background. The likelihood function has the form 



N 



C \ \ [f s - P a (mi,ti,<Tu)- 



i=i 



(l-f.)-Pd,(mi,ti,<rti)]. (11) 



The parameter f s denotes the fraction of signal B® — > 
J/tpfo decays and P s and P c b the probability density 
function of signal and combinatorial background, respec- 
tively. To enhance the signal-to-background ratio in the 
selected sample, we use only B° candidates with decay 
times larger than 0.2 mm/c = 0.67 ps. This requirement 
suppresses background by a factor of 40 and reduces the 
prompt background component to a negligible level while 
keeping about two thirds of the signal events. 

The B® signal mass PDF is parametrized as for the 
branching ratio measurement. The PDF in decay time 
is parametrized with an exponential function convolved 



with a Gaussian resolution function. The width of the 
Gaussian is given by the candidate-specific estimated de- 
cay time uncertainty au scaled by a common factor St 
which accounts for possible discrepancies between esti- 
mated and actual resolutions. The scaling factor St is 
determined in a fit to data dominated by prompt back- 
ground, selected by requiring < t < 0.3 ps. In the final 
fit, St is a free parameter with a Gaussian constraint in- 
cluded as additional factor in the likelihood in Eq. (jlll) . 
The PDF in decay time uncertainty is parametrized by 
an empirical function. We use a log-normal distribution 
with parameters [i, 9, and a defined as 



T>((Xti\fx,0,a) 



1 



2na(ati - /x) 



(12) 



for au > [i and zero otherwise. Given the rather small 
statistics of the B® signal we derive the parameters using 
simulated B° s — > J/ipfo events and Gaussian constrain the 
values in the fit to data. The widths of the Gaussian con- 
straints are chosen to cover possible differences between 
simulation and data. 

The combinatorial background is described by two 
components, a long-lived part for the background from 
6-hadron decays and a short-lived part for the tail from 
mis-reconstructed prompt events. The mass PDF is com- 
mon to both components and parametrized by a linear 
function. The decay time PDF of each component is de- 
scribed by an exponential convolved with the same reso- 
lution function as used for signal. Both decay time uncer- 
tainty PDFs are again modeled using log-normal distri- 
butions. The parameters of each log normal distribution 
are independent of the distribution of the B® signal. 

All parameters of the combinatorial background are 
determined from the fit. The yield, the mass resolution 
scale factor, and the lifetime of the B® signal are also 
left to float freely. The decay time uncertainty parame- 
ters of the signal and the resolution scale parameter are 
Gaussian constrained. Using an ensemble of simulated 
experiments we verify within 1% that the fit is unbiased 
and returns proper uncertainties. 



B. Systematic uncertainties 

We investigate several possible sources of systematic 
uncertainties. These are broadly separable into two 
classes: the first dealing with the parametrization of the 
PDFs and the second with possible biases in the selection 
or reconstruction. 

We first investigate our assumption of the mass shape 
of combinatorial background. We determine the relative 
change of the B® lifetime between a fit with a first and 
a third order polynomial background mass model. For 
fits in different invariant mass ranges, we find an average 
difference of 0.010 ps, which we assign as the system- 
atic uncertainty. The systematic uncertainty assigned to 
the signal mass shape has contributions from the limited 
knowledge of the mean position and from the assumed 
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shape parametrization. Both effects are evafuated in the 
same way as for the branching ratio measurement and 
yield a systematic uncertainty of 0.005 ps. There are 
two assumptions made for the decay time PDFs; one is 
the resolution scale factor, St, which is known only with 
limited precision and the other is the shape of the com- 
binatorial background. The uncertainty of the scale St 
is included directly in the statistical uncertainty of the 
fit as the parameter is allowed to vary within a Gaussian 
constraint. To quantify the size of the contribution, we 
repeat the fit with St fixed to its central value and find 
the quadratic difference in uncertainty to the original fit 
to be 0.005 ps. To estimate the effect of the assumed de- 
cay time PDF of combinatorial background, we employ 
an alternative fit method which does not need a decay 
time parametrization of the background. We split the 
data into 20 decay time bins and simultaneously fit the 
invariant mass distributions with independent parame- 
ters for the background in each bin. The signal yield per 
bin is given by the total signal yield times the integral of 
the signal decay time PDF over the time bin, where the 
same PDF parametrization as in the default fit is used. 
The difference in the fit results is taken as a measure 
of the systematic uncertainty due to the background de- 
cay time PDF. To avoid possible statistical fluctuations 
in this estimate we repeat the comparison for different 
selection requirements and assign the average difference 
of 0.021 ps as systematic uncertainty. The third kind 
of systematic effect addresses the uncertainty of the at 
PDFs. The main effect is the distribution for signal de- 
rived from simulated events. The uncertainty is already 
included in the statistical error since the parameters are 
Gaussian constrained in the fit. The contribution due 
to modeling of the decay time uncertainty distribution, 
estimated from a comparison of fit results with fixed and 
constrained parameters, is 0.015 ps. 



TABLE II: Summary of assigned systematic uncertainties for 
the lifetime measurement. The uncertainties in parentheses 
are included in the statistical uncertainty via Gaussian con- 
straints in the fit. 



Source 


Uncertainty [ps] 


Background mass model 


0.010 


Signal mass model 


0.005 


Decay time uncertainty scale 


(0.005) 


Background decay time model 


0.021 


Decay time uncertainty model 


(0.015) 


SVX alignment 


0.007 


Total 


0.03 
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FIG. 9: (color online) Invariant mass distribution with fit 
projection overlaid. 



For the second class, we verify that our candidate se- 
lection does not introduce any significant bias. A bias 
in the mass distribution could artificially enhance or de- 
crease the amount of signal candidates while a bias in 
decay time could directly affect the extracted lifetime. 
We verify on a background-enriched sample selected by 
requiring t < 0.01 cm/c that no artificial peak is ob- 
served for any neural network output requirement. With 
a high statistics sample of simulated events we check that 
the selection does not bias the fitted lifetime. A possible 
lifetime bias introduced by the trigger has been stud- 
ied in a previous CDF analysis 45[ and is negligible in 
our measurement. Finally the alignment of the tracking 
detectors is known only with finite precision. Previous 
measurements found that the uncertainty on the lifetime 
due to a possible misalignment is 0.007 ps [45 1. 



All the contributions are added in quadrature and yield 
a total systematic error on the lifetime of 0.03 ps (1.5%). 
A summary of the systematic uncertainties on the life- 
time is provided in Table [TTJ 



C. Lifetime result 

Performing the likelihood fit to the selected data we 
extract the B® lifetime in B® —> J/ifjfo decays 

t{B° s Jyty/o) = 1.70±g:J?(stat) ± 0.03(syst) ps. (13) 

In Figs. [9] to [11] we show the data together with the pro- 
jection of the fit. 



VI. CONCLUSIONS 

We confirm the observation of the B® — > J/0 /o (980) 
decay from the LHCb [2l| and Belle [22[ experiments. 
The observed signal is the world's largest and we perform 
the most precise measurement of the ratio of branch- 
ing fractions Rf a /^ between B® —> J/ipfo an d B° — ► 
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FIG. 10: (color online) Decay time distribution with fit pro- 
jection overlaid. 
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FIG. 11: (color online) Decay time uncertainty distribution 
with fit projection overlaid. 



J/tpcj)(980) decays: 

= g(gg -> J/^/o(980)) g(/ (980) -> TT+vr-) _ 

0.257 ± 0.020(stat) ± 0.014(syst). (14) 

In this result we assume that the observed signal is solely 
due to the decay B° s — > J/ip fo (980) and correct for the 
acceptance of the invariant mass selection of the pion 
pair. Using the world average B® — > J/tp<p branching 
fraction [lj] Rf /,p can be converted into the product of 
branching fractions of 

B(B° S -> J/^/ (980))6(/o(980) tt+tt-) = 

(1.63 ± 0.12 ± 0.09 ± 0.50) x 10" 4 , (15) 



where the first uncertainty is statistical, the second is 
systematic, and the third one is due to the uncertainty 
on the B® — > J/4>4> an d 4> — > K + K~ branching fractions. 
The measurement presented here agrees well with the 
previous measurements of this quantity and with theo- 
retical predictions. 

Moreover, our sample allows us to measure the B® life- 
time in the B® — > J/ip fo (980) decay mode: 

t{B° ->• J/0/ o (98O)) = 1.70t^i(stat) ± 0.03(syst) ps. 

This is the first measurement of the B® lifetime in a decay 
to a pure CP eigenstate. In the context of the standard 
model the lifetime measured in this decay mode to a CP- 
odd final state can be interpreted as the lifetime of the 
heavy B® eigenstate. The measured value agrees well 
both with the standard model expectation as well as with 
other experimental determinations. 

While the precision of the lifetime measurement is 
still limited by statistics, it provides an important cross- 
check on the result determined in B® — > J/tp4> decays, 
which relies on an angular separation of two CP eigen- 
states. Furthermore, the measured lifetime can be used 
as an external constraint in the B® — » J/ip4> analysis 
to improve the determination of the CP- violating phase 
in the B® — > J/ifxj) decay. The lifetime measurement 
in B® — > J/!/)/o(980) decays is also the next step to- 
wards a tagged time dependent CP-violation measure- 
ment, which can provide an independent constraint on 
the CP violation in B® mixing. 
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